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ABSTRACT 

Maunder Minimum forms an archetype for the Grand minima, and detailed knowledge of its 
temporal development has important consequences for the solar dynamo theory dealing with long- 
term solar activity evolution. Here we reconsider the current paradigm of the Grand minimum 
general scenario by using newly recovered sunspot observations by G. Marcgraf and revising 
some earlier uncertain data for the period 1636-1642, i.e., one solar cycle before the beginning 
of the Maunder Minimum. The new and revised data dramatically change the magnitude of 
the sunspot cycle just before the Maunder Minimum, from 60-70 down to about 20, implying a 
possibly gradual onset of the Minimum with reduced activity started two cycles before it. This 
revised scenario of the Maunder Minimum changes, through the paradigm for Grand solar /stellar 
activity minima, the observational constraint on the solar/stellar dynamo theories focused on 
long-term studies and occurrence of Grand minima. 

Subject headings: Sun; activity — sunspots — dynamo 

1. Introduction lar activity, which is grossly dominated by the 

quasi-periodic wll-year cycle. On top of that. 
The Sun is the only star whose features can be secular variability is superposed ranging from very 

studied in great detail and on long time scale, ^igh activity during the late 20-th centu ry down 



forming a paradigm for a large population of ^o ver y quiet per iods of Grand minima (lUsoskin 



sun-like' stars. The Sun exhibits a great deal [200! iHathawavl [2010h . Generally, the 11-year 



of magnetic variability generically called the so- cyclicity is understood in terms of solar dynamo 
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( Charbonneau 201Ci) . while the secular/mihenial 
variabihty still remains unclear. A particular 
enigma for solar/stellar dynamo is the occurrence 
of Grand minima, that is not an intrinsic feature of 



the standard mean- field dynamo (jSokolofn 120041) . 
Accordingly, the occurrence of Grand minima is 
often modelled by an ad hoc approach, including, 
e.g., stochastic o r chaotic driven processes (e.g. 
Moss et aD l2008l ) . Therefore, observational con- 
straints on solar/stellar dynamo with respect to 
Grand minima are crucially important. While the 
statistic of Grand minima/maxima of solar activ- 
ity is studied t o some extent using the cosmogenic 
isotope data (lEddyl Il977l : IStuiver fc Braziunas 



1989t lUsoskin et al-lboOTt lAbreu et al.ll2008l ). the 

variability of solar activity during and around a 
Grand minimum is not precisely known. This is 
caused by the fact that presently we have more 
or less detailed information on only one Grand 
minimum - the Mau nder rninimu m in the second 
half of 17-th century (|Eddvll 19761) . which serves as 
an archetype for Grand minima in general. The 
present parad igm for a general scenario of a Grand 
mmimum fsee lUsoskin fc Mursulall2003l and refer- 
ences therein) is widely considered as a constraint 
for solar dynamo models. Therefore, every piece 
of information on solar activity during that period 
is extremely valuable. 

Temporal variability of solar activity around 
the Maunder minimum is usually studied using 
histor ical telescopic observations of sunspots since 
1610 (|Vaauero fc Vazauez|[2009l) . The period be- 
fore the Maunder minimum (i.e., the first half of 
the 17-th century) was very uncertain in the origi- 
nal Wolf sunspot series, but the data were greatly 



i mpro ved after tremendous work of iHovt et al 



(|l994l ) and lHovt fc SchattenI (|l998l ). who collected 
large amount of historical sunspot records, form- 
ing the group sunspot number (GSN) series. The 
solar activity variability before 1650 AD is shown 
in Fig. [T] according to the official GSN series 
(dotted curve) and r eanalyzed with impr oved sta- 
tistical techniques ( Usoskin et al.l l2003l - grey 
curve). The first observed solar cycle 1610-1618 
was pretty high (maximum annual GSN values 
above 100), and this is quite reliable as it was cov- 
ered by the direct data inclu ding sunspot drawings 
and counts (jHovt &: Schattc n 199§). The second 
cycle 1618-1630 was lower (maximum GSN values 
30-40), but the quality of data is not good. The 



next cycle 1635-1645 is marked as high in the GSN 
series, and this gave rise to the idea of an abrupt 
onset of the general Grand minimum scenario. 

However, the GSN series is also not complete, 
and new data, which remained unnoticed by Wolf 
and his successors including Hoyt & Schatten, 
are continuously recovered in v arious places, of- 
ten outside major observato ries ('Arlt"2008 [ 120091: 



Casa.s et "all I2OO6I : IVaauero 2004; .Vaauero et al 
2OO5I l2007al ). This leads to a revision of some 



parts of the GSN series and requires re-analysis of 
some results. Here we reconsider the paradigm of 
Grand minimum general scenario by using newly 
recovered sunspot observations as well as revising 
some earlier uncertain data for the period 1636- 
1642, i.e., one solar cycle before the beginning of 
the Maunder minimum. 

2. Revised sunspot data for 1636 1642 



As discussed in great detail bv lVaauerol(l2007bh 
sunspot data in 1610-1644 (before the Maun- 
der Minimum) contains numerous gaps and un- 
certainties. This problem is specially acute for 
the period from 1636 through 1642. For that 
period the following data sets were available in 
the GSN series: (1) three years (1636, 1637 and 
1641) without any solar observational record; (2) 
two years (1640 and 1642) with little amount of 
records by Horrox, Gassendi, Hevelius, Scheiner 
and Rheita (2 and 37, respectively); and (3) two 
years (1638 and 1639) with estimated data based 
on Crabtree data. Here we use the GSN database 
(ftp:/ /ftp.ngdc.noaa.gov/STP /SOL ARJDATA / 
SUNSPOT.NUMBERS/GROUP_SUNSPOT.NUMBERS/), 
but critically revise it for 1636-1642. 

First we add newly recovered sunspot data by 
G. Marcgraf which are not included into the GSN 
database. Georg Marcgraf (1610 -1648) was a Ger- 
man naturalist and astronomer (jN^ligii). We 
have consulted his manuscripts of astronomical ob- 
servations preserved in Leiden Regional Archive 
(Collection Marcgraf LB7000/1) and in the Na- 
tional Library of Portugal (Mss. 6, n° 37). We 
have recovered sunspot records for the year 1637, 
just the year with no records in the GSN database. 
These records are preserved in the 'Collection 
Marcgraf (Leiden Regional Archive, LB7000/1) 
in the documents labelled as 22, 24a, 24b and 52 
according to the numeration made by the historian 
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of science John North in the 1980s (North 1989). 
Fig.[2]depicts an example of solar disc drawings by 
Marcgraf showing one sunspot from 9 to 12 June 
1637 (document 24a). Sunspot records by Marc- 
graf cover 17 days of 1637, that were included in 
the new database (see Table [T]). 

Another addition is correction of Horrox's 
records. Jeremiah Horrocks or Horrox (1618- 
1641) was an English astronomer who predicted, 
observed and recorded the transit of Venus in 
1639. His sunspot records appeared in his 'Opera 
Posthuma' (HorrocciiJ673 ). Horrox used the Ju- 
lian Calendar and, therefore, we have c onverted 
these dates to Gregorian Calendar ( Vaquerol 
2007bl ). Moreover, Horrox noted "Maculae duas 
in Sole" (two spots in the Sun) in 30 Oct - 1 
Nov 1638. We have interpreted these two spots as 
one sunspot group to make this record compatible 
with the only group observed by Gassendi on the 
same dates. 

Next we noted that the filling of the years 1638 
and 1639 in GSN was done base d on an unclear es- 
timate rather than on real data. iHovt &: Schatteii 
(|l998l ) wrote in their Bibliography appendix: " Ac- 
cording to a letter by Crabtree the average number 
of spot groups seen in 1638 and 1639 were 4-5 per 
day. The database has Greenwich fill values to 
give 4-5 groups per day. This substitution tech- 
nique was used to simplify the analysis. This is 
the only place in the entire database where we do 
this type of substitution." However, the number of 
the actually observed groups was smaller than the 
number of estimated groups for these years in th e 
reconstruction (see Figure 5 in Vaquerol 2007b). 
Therefore, we decided not to use this non-observed 
values for our estimations and eliminate Crabtree 
data. 

We have also made some minor corrections to 
the database using the original sources. We have 
eliminated one spurious observation by Gassendi 
on 1 Dec 1638, because this record doe s not ap- 
pear in his astronomical observations ([Gassendi 
16581 ). Moreover, we have incorporated one 
sunspot record by Horrox in 4 December 1639. 
Curiously, the treatise "Venus in sole visa" [Venus 
in transit across the Sun] by Horrox was published 
by Johannes H evelius at his own expense in 1662 
( Heveliu^ 1662 ). In this work, Horrox noted that 
he saw one small and common spot i n the Sun on 
4 Dec 1639 (p. 115 of lHeveliusl[l662[) . Finally, we 



have also changed the record by Rheita in 1642 
using the original source (pp. 242-243 of [Rheita 
1645li . 

All these additions and changes in the original 
Hoyt and Schatten database are listed in Table [1] 
and summarized below: 

1. Newly recovered sunspot records by G. Mar- 
cgraf are added; 

2. The estimated (not observed) values from 
Crabtree's comments (1638-1639) are dis- 
carded; 

3. Dates and numbers of sunspot groups are 
corrected for Horrox's observations; 

4. A spurious observation by Gassendi on 1 De- 
cember 1638 is eliminated; 

5. One missing sunspot record by Horrox in 4 
December 1639 is added. 

6. The record by Rheita in 1642 is corrected. 

Using the revised data presented in Table [U we 
have evaluated the annual sunspot numbers for 
the period of 1637-16 42, emp lo ying t he sta tistical 
method proposed by lUsoskin etaU (|2003l) . The 
method is based on comparing the actual sparse 
data (sample population) to the daily sunspot 
data in 1850-1996 (reference population), as- 
suming constancy of the statistical properties of 
sunspot activity. For a given sample population 
of daily measurements within a month, months 
in the reference population are found that con- 
tain the same subset of daily values. A statistical 
distribution of the corresponding monthly means 
is then built that allows to estimate the mean 
and uncertainty of monthly sunspot numbers, re- 
constructed from sparse daily o bservations (see 
full details in lUsoskin et al.ll2003[ ). From monthly 
mean values, yearly values can be obtained in the 
same way. The newly computed yearly sunspot 
numbers for 1637-1642 are given in Table [2] and 
shown in Fig. [T] as thick solid line with error bars. 

As one can see from Fig. [U the revised and up- 
dated data for 1637-1642 has essentially changed 
the profile of temporal variability of sunspot data 
before the Maunder minimum. In particular, the 
last solar cycle before the minimum now appears 
quite modest, with the peak value of 20 ± 15 com- 
pared to about 70 in the GSN series. The new 
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scenario, in accordance with the new database, im- 
plies low solar activity roughly two solar cycles be- 
fore the beginning of the Maunder minimum. This 
suggests that transition from the normal activity 
to the deep mi nimum was no t as sudden as previ- 
ously thought (IUsoskinll2008l ). and the descent of 
solar activity might have started as early as in the 
1610s. 

3. Discussion and Conclusions 

Since the Maunder minimum forms an archetype 
for the Grand minima, detailed knowledge of 
its temporal development has important conse- 
quences for the solar dynamo theory dealing with 
long-term solar activity evolution. The ge neral 
dynamo theory (see, e.g. Charbonneau 201Clt ) can- 



not naturally reproduce occurrence of Grand min- 
ima and requires some prescribed changes in the 



the Maunder minimum ( 


3.g. Vitinskv et al. 


1986 


Sokoloff & Nesme-Ribes' 


1994; 'Frick et al.l 


1997 


Usoskin et al...2000,; Usoskin.,2008) is that transi- 



tion from the normal high activity to the deep min- 
imum was sudden (within a few years) and with- 
out any apparent precursor, while the recovery to 
the normal activity level was gradual, taking sev- 
eral decades. The abrupt onset of a Grand minima 
forms a strong constraint, as only few models with 
stochastically driving fo rces can 'naturally' pro- 
duce such a feature (e.g. ICharbonneau fc Dikpati 
[2000), while others require special ad hoc assump- 
tions. Presently, several models can reproduce, 
with different approaches, the proposed scenario of 
a Gr a nd minimum (e.g. [ Charbonneau &: Dikpati 



2000: ICharbonneau et 



aLri2004l: 



2009at karaklbOlOl: IPassos fc Lopesll201ll) . 



Usoskin et al 



On the other hand, many dynamo models suc- 
ceed in predicting Grand minima to start gradu- 
ally, through a continuous decrease of the activity 
level to the deep minimum, followe d by a grad- 



ual re c overy. Just to rnention a few. | Kiiker et al 



(11999') : 'Weiss fc Tobia s (2000*1: 'Brooke et al.l (l2002l) : 
feushby (20061 ) : lMoss"et al., (gOOi). Thus, it is cru- 
cial to know the solar activity evolution before the 
Maunder Minimum with high confidence. 

The revised sunspot records presented here im- 
ply that the scenario of the Maunder Minimum 
can correspond to a gradual onset, contrary to the 
earlier consideration, thus affecting observational 



constraints on the solar/stellar dynamo. Unfortu- 
nately, cosmogeni c isotope data can hardly resolve 
individual cycles ( Steig et "al]ll996t lUsoskin et al 



r2009bl) to clearly answer this question, but a cau- 
tious statistical study of radiocarbon ^''^C data in 
tree rings for the Maunder and Sporer solar min- 
ima indicate a possible lengthening and attenua- 
tion of solar cycles a f ew decades b efore the onset 
of a Grand minimum |Miyahara et al. 2008). 

The major findings of this work can be summa- 
rized as follows. 

1. Using newly recovered sunspot records by 
Georg Marcgraf and carefully revised data 
for other observations (Table |T]) , we provide 
a new sunspot number series (Table [2]) for 
the period 1636-1642. 

2. The new data dramatically change the mag- 
nitude of the sunspot cycle just before the 
Maunder Minimum, from 60-70 down to 
about 20, implying a possibly gradual on- 
set of the Minimum with reduced activity 
started two cycles before it. 

3. This revised scenario of the Maunder Min- 
imum changes, through the paradigm for 
Grand solar/stellar activity minima, the ob- 
servational constraint on the solar/stellar 
dynamo theories focused on long-term stud- 
ies and occurrence of Grand minima. 

Thus, we have essentially revised the sunspot 
data prior to the Maunder Minimum leading to the 
revisited observational scenario of a Grand mini- 
mum of solar activity. The present results are ex- 
pected to impact development of models dealing 
with long-term solar/stellar activity evolution. 
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Support from the Junta de Extremadura and Min- 
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from the Finnish Academy. 
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Table 1: Sunspot observations (date, number of 
sunspot groups G and the observer) for 1636-1642 
used in this stay. The last column comments 
on the relation of the observation to the HS98 



(iHovt fc Schatten ''l998'l database. 



Date 


G 


Observers 


Comment 


19-20 Jan 1637 


1 


Marcgraf 


New data 


22 Jan 1637 





Marcgraf 


New data 


5 Feb 1637 


1 


Marcgraf 


New data 


9-12 Jun 1637 


1 


Marcgraf 


New data 


13 Jun 1637 





Marcgraf 


New data 


21 Sep 1637 


2 


Marcgraf 


New data 


24-25 Sep 1637 


2 


Marcgraf 


New data 


28 Sep 1637 





Marcgraf 


New data 


10 Oct 1637 





Marcgraf 


New data 


12-13 Oct 1637 


1 


Marcgraf 


New data 


15 Oct 1637 





Marcgraf 


New data 


1-3 Jun 1638 


2 


Horrox 


Corrected date 


29 Octl638 





Gassendi 


HS98 


30 Oct-1 Nov 1638 


1 


Horrox 


Corrected date and G 


30 Oct-1 Nov 1638 


1 


Gassendi 


HS98 


4 Dec 1639 


1 


Horrox 


New data 


21-22 Aug 1640 


1 


Scheiner 


HS98 


9-22 Jun 1642 


1 


Rheita 


Corrected date and G 


26 Oct 1642 





Hevelius 


HS98 


28 Oct 1642 


1 


Hevelius 


HS98 


31 Oct-1 Nov 1642 


1 


Hevelius 


HS98 


3-4 Nov 1642 


1 


Hevelius 


HS98 


6 Nov 1642 


3 


Hevelius 


HS98 


8 Nov 1642 


2 


Hevelius 


HS98 


9 and 11-17 Nov 1642 


1 


Hevelius 


HS98 


18 Nov 1642 





Hevelius 


HS98 



Table 2: Yearly group su nspot numbers for 1637- 
1642: formal values Rg (|Hovt fc SchattenI Il998l) 
as well as the calculated here weighted values Rw 
with ±g uncertainties. 
Year R„ Rw 



1637 


n.a. 


13±2 


1638 


68.7 


19±6 


1639 


76.8 


21±15 


1640 


15 


17±12 


1641 


n.a. 


n.a. 


1642 


47.3 


13±3 
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140 



120 



1610 



- Hoyt & Schatten (1998) 
-Usoskinetal. (2003) 
-This work 




1620 



1630 



1640 Year 



Fig. 1. — Annual sunspot numbers in the first 
half of 17-th century. G roup sunspot numbers Rg 
(jHovt &: Schattenlll998| ) are shown by the dotted 
Une, weig hted sunspot number , based on the same 
data set ( Usoskin et al. 20031 ) by grey line, and 
the weighted sunspot number estimated in this 
work by black line. 




Fig. 2. — An example of the solar disc drawing 
by Marcgraft depicting one sunspot from 9 to 12 
June 1637. 
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